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Many electronic systems (e.g., the cuprate superconductors and heavy fermions) exhibit striking 
features in their dynamical response over a prominent range of experimental parameters. While there 
are some empirical suggestions of particular increasing length scales that accompany such transi- 
tions in some cases, this identification is not universal and in numerous instances no large correlation 
length is evident. Using a correspondence between dissipative classical and quantum many-body 
systems, we show that, even in the absence of imposed disorder, many continuum systems (and 
possible lattice counterparts) may exhibit zero-point "quantum dynamical heterogeneities" wherein 
the dynamics, at a given instant, is spatially non-uniform. Towards this end, we extend a known 
mapping between finite temperature classical Fokker-Planck systems and quantum systems at zero 
temperature to include general non-equilibrium dynamics. While the static length scales accompa- 
nying this phenomenon do not seem to exhibit a clear divergence in standard correlation functions, 
the length scale of the dynamical heterogeneities can increase dramatically. We furthermore show 
how a hard core bosonic system can undergo a zero temperature quantum critical metal-to-insulator- 
type transition with an extremely large effective dynamical exponent z > 4 that is consistent with 
length scales that increase far more slowly than the relaxation time as a putative critical transition 
is approached. We suggest ways to analyze experimental data in order to adduce such phenomena. 
Our approach may be applied to quenched quantum systems. 

PACS numbers: 05.30.-d, 03.67.Pp, 05.30. Pr, 11.15.-q 



I. INTRODUCTION 



A prominent centerpiece in the understanding of nu- 
merous systems is Landau Fermi-liquid theory (LFL); 
this theory allows the understanding of phenomena such 
as conventional metals and low temperature He-3 liq- 
uids. LFL is centered on the premise that the low en- 
ergy states of interacting electron systems may be cap- 
tured by long-lived fermionic quasi-particles with renor- 
malized parameters (e.g., effective masses that differ from 
those of the bare electron). The last three decades have 
seen the discovery of materials in which electronic be- 
havior deviates from simple LFL. These "singular" or 
"non-Fermi liquids" (NFL) include the pseudo-gap re- 
gion of the high-temperature cuprate superconductors 
and "heavy fermions" (in which, as befits their name, the 
effective electron mass becomes very large). While there 
are clear indications of changes in the dynamics in t hese 
systems-including putative quantum critical points^L 
there is, in most cases, no clear experimentally measured 
length scale that exhibits a clear divergence. A quantum 
critical point is associated with a continuous phase transi- 
tion at (absolute) zero temperature. Typically, this may 
occur in a system whose transition temperature is driven 
to zero by doping or the application of magnetic fields 
or pressure. Within a quantum critical regime, response 
functions follow universal power law scaling in both space 
and time. Specifically, at a quantum critical point, the 
effective infrared (IR) fixed point theory exhibits scaling 



invariance in space-time: t —¥ Xt, x — > X x ^ z x with a dy- 
namical exponent z that can, depending on the theory at 
hand, assume various canonical values. Unlike classical 
critical points whose associated critical fluctuations are 
confined to a narrow region near the phase transition, 
quantum critical fluctuations appear over a wide range 
of temperatures above the quantum critical point. These 
fluctuations may generally lead to a radical departure 
of the system's electronic properties from standard LFL 
type behavior. These features are anticipated to be com- 
mon across many strongly correlated electronic systems 
and may be associated, in some electronic systems, with 
a change of Fermi surface topology®. 

Many NFL systems exhibit numerous phases (includ- 
ing, quite notably, superconductivity). Indeed, compet- 
ing orders and proliferation of multiple low energy states 
can lead to spatially non-uniform glassy characteristics® 
and associated first order quantum transitions®. The 
length scales characterizing these electronic systems un- 
dergo a much milder change than the corresponding 
changes in the dynamics. All of this suggests an effec- 
tive infrared (IR) fixed point quantum field theory is in- 
variant under scaling in time but not in space- i.e., the 
effective dynamical exponent z — > oo. Response func- 
tions such as those of the marginal Fermi liquid form 
describing cuprates near optimal doping show a marked 
frequency dependence but essentially no spatial momen- 
tum dependence. In this work, we will suggest how many 
quantum systems might exhibit very large effective dy- 
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namical exponents. 

It is natural to first look elsewhere in physics where 
similar phenomena appear. One arena immediately 
comes to mind. In classical structural glasses there is 
a dramatic change in the dynamics as a supercooled liq- 
uid is quenched into a glass without the appearance of 
easily discernible large changes in measurable standard 
static length scales. While the ergodicity breaking that 
accompanies a glass transition cannot occur in a finite 
size system, it essentially mandates the appearance of a 
diverging static length scale^, but such divergent length 
scales generally do not simply manifest themselves in 
bare standard correlation functions. General correlation 
functions which may monitor subtle changes include the 
"point to set"^ correlations and others. Practically, in 
most instances^ no clear signatures of divergent length 
scales are easily seen in standard static two-point corre- 
lation functions. 

A far more transparent growth in length scales is seen 
from four-point correlation functions that quantify the 
change in correlations as the system evolves in time. 
These correlation functions afford a glimpse into the 
length and time scaling which describe dynamical het- 
erogeneities that characterize the spatially non-uniform 
rate of change or dynamics in the system. The length 
scale associated with these heterogeneities was seen to 
grow as the characteristic relaxation times increased. 

We may use similar correlation functions to character- 
ize strongly correlated electronic systems in which there 
are strongly discernible changes in the dynamics but no 
obvious experimentally accessible tools that point to ac- 
companying divergent length scales in a general waj^. 
To our knowledge, to date, dynamical heterogeneities 
(nor general static measures such as those of the point- 
to-set method) have not been systematically probed for 
in electronic systems nor has their existence been estab- 
lished as a matter of principle in quantum system. Ini- 
tial ideas concerning non-uniform doping-driven hetero- 
geneities were discussed irP^. In this work, we flesh out 
the blueprint for a proof outlined, by one of us, irP^ and 
we will provide concrete "matter of principle" theoret- 
ical testimony to the emergence of quantum dynamical 
heterogeneities in clean systems and related properties in 
quantum many body systems. Some time after^, the au- 
thors of^ have further confirmed the existence of quan- 
tum dynamical heterogeneities in certain dissipative spin 
systems. In what follows, we focus on systems such as 
structural glasses or jammed systems with no external 
disorder. 

As is well appreciated, such translationally invariant 
systems that are free of external disorder may, on their 
own, display non-uniform spatial patterns concomitant 
with interesting dynamical properties. For times far 
shorter than the equilibration times, the two-point auto- 
correlation function records dynamic fluctuations at all 
points x in d dimensional space, 



forms an analogue of the Edwards- Anderson^ order pa- 
rameter that appears in spin glasses. The spatial corre- 
lation amongst pair products of time separated products 
of fields (such as those in Eq. ([I])) at different spatial 
sites is a four-point correlation functional that attempts 
to measure cooperation 

Giix-y.t) = {5cp{x,t)54>{x,Q)54>{y,t)84>(y,Q)) - C 2 (t). 

This four-point correlation function relates the dynamics 
at two different spatial points x and y; due to transla- 
tional invariance, this correlation function depends only 
on (x — y) and not on x and y separately. 



II. REVIEW OF A RELATION BETWEEN 
VISCOUS CLASSICAL AND QUANTUM MANY 
BODY SYSTEMS AND A NEW 
GENERALIZATION 

In order to illustrate how, as a matter of principle, 
the physics of such classical dissipative systems can ap- 
pear in clean quantum system s at zer o temperature, we 
employ and extend a mappin g 15 * 16 ! 20 " 22 ! between classi- 
cal dissipative systems and quantum many body systems 
to include general dynamical (including non-equilibrium) 
systems. Related work concerning dynamics in Rokshar- 
KivelsorP^l systems also appears irPl Below, we briefly 
review this mapping. We will then discuss some of its 
(new) consequences. Details of the derivation of our new 
result concerning a classical to quantum correspondence 
for the dynamics for general dissipative classical systems 
have been relegated to Appendix [C] 

The crux of the mapping between dis sipative classical 
systems and many body bosonic theorie d 15 * 16 * 20 ! is the re- 
alization that the equation of motion for a dissipative 
(or "Aristotelian" ) classical system is a first order differ- 
ential equation in time just as the Schrodinger equation 
is. Using this equivalence, systems obeying the Langevin 
equation, 



7i-Jr = -ViVn(xi, 



;X N ) + ffi(t), 



(3) 



[with i the particle index, 7^ the coefficients of friction, 
rjf (t) the Gaussian noise with 



(V?(t)r]-(t')) = 2T clli 8 ij 5 a p8{t - t'), 



(4) 



where T c ; is the effective temperature of the classi- 
cal sys tem, an d, a, (3 — 1,2, ...,c£] can be exactly 
mapped 15 * 16 l 2 ° l onto a quantum many body system of 
bosons with effective mass m, = r ) i /{2T c i) at zero tem- 
perature which is governed by the Hamiltonian 



H 



C(t) = / d d x(5(j)(x,t)8(l>(x,0)), 



(1) 
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We will term the quantum system of Eq. |5]) as the 
"dual quantum system" associated with the classical sys- 
tem of Eqs. ([3| [4]). The many body quantum poten- 
tial VQuantum({^}) is constructed from the gradients of 
the classical potential energy Vn as in Eq. Un- 
der this mappingp^HIH a dissipative classical system with 
a potential energy Vn that captures repulsive hard core 
spheres maps onto a quantum system at zero temper- 
ature with (as is apparent in the many body potential 
energy VQ ua ntum) similar dominant hard-core interactions 
(augmen ted by s oft sticky interactions)!^. Earlier work 
advance d 15 * 17 * ^ the rudiments of this mapping and fur- 
ther suggested dynamical aspects that might follow from 
it. In the current work, we explicitly derive and prove 
a general and rather powerful relation between arbitrary 
correlation functions in dissipative classical systems with 
time varying potentials (necessary for our discussion of 
quenching) and relate these to corresponding correlation 
functions in zero-temperature quantum systems. In the 
below we provide a key result of this m apping. In Ap- 
pendix [cj we summarize, followin ^ 15 * 17 * ^, more detailed 
aspects of this mapping. We now proceed to set the stage 
for our new result and its consequences. [We should fur- 
ther note that although we will focus on the mapping 
from classical systems to corresponding quantum ones, it 
is also possible to go in the opposite direction and map 
quantum mechanical systems with known non-degenerate 
ground states onto classical dissipative systems (see ap- 
pendix [H]).] In what follows, we consider a general clas- 
sical two-point correlation function of the form 



G 



classical 



(t) = (O(t)O(0)), 



(6) 



where 0(t) is any quantity at times t > 0. A notable 
extension of our work, presented in Appendix [C] is that 
our mapping covers general non-equilibrium time depen- 
dent Hamiltonians in which only the initial (or the final) 
classical state is in thermal equilibrium at an initial (or fi- 
nal) temperature. In such a case, the corresponding dual 
quantum problem evolves unitarily with a time depen- 
dent Hamiltonian. Specifically, as we will elaborate on 
in Appendix [Cj we find a very simple and general result 
concerning a rather simple "complexification" of the time 
coordinate (or Wick-type rotation) in going from general 
classical correlation functions of the form of Eq. ^ to 
their corresponding bosonic counterparts. The final re- 
sult is rather simple (Eq.Q). For any quantity 0(t) if the 
quantum correlation function GQ uantum ((f, t) is evaluated 
with th e qua ntum many body potential VQ uan tum({^}) 
of Eq. ( A12 1 while the corresponding classical correla- 



tion function is computed for a system with a potential 
Classical ({%}) then a very simple relation exists, 



-RQuantum(i) = (Classical - Cclassical (**)) ■ ( 7 ) 



Eq. for the response function provides a new entry 
in the mapping between the finite temperature classi- 
cal system of Eq. ^ governed by a potential Vn and 
its corresponding quantum zero-temperature dual with a 



quantum many body potential energy VQuantum ({a?}) in 
the Hamiltonian of Eq. ([5| . 

Eq. is a key result of this work. In Appendix [C] we 
provide a detailed derivation of Eq. |7]). More general 
relation for correlation functions in systems with time 
dependent Hamiltonians (concerning evolving the system 
out of equilibrium classically or out of its ground state 



quantum mechanically) appear in Eqs. (C16 C17 C18|. 
Typically, in glassy systems, the correlation function of 
Eq. (JsJ) is a superposition of many decaying modes- This 
distribution of modes will manifest as a distribution of os- 
cillatory modes in the corresponding dual quantum prob- 
lem. In many cases, this will lead to zero temperature 
quantum dynamics of the dual system that, with addi- 
tional oscillations, will emulate the finite classical dy- 
namics. For instance, if, as it precisely occurs in viscous 
systems with overdamped dynamics (for which Eq. ^ 
applies), if for times t > 0, 



G c | assica i(i) =^exp[-(t/T) Q ] 



(8) 



then, correspondingly for all of these positive times t (see 
Appendix |e| , 
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(9) 



With the aid of the general relation of Eq. ([7| , the quan- 
tum correlation function that corresponds to a general 
stretched exponential correlation function in the classical 
arena can be computed analytically and is, indeed, given 
by Eq.(|9|. In general, there are, of course, subtle issues 
relating to the analytic continuation of general functions 
to imaginary time. A trivial yet important particular re- 
alization of Eq. ([7| is that static correlations (t = 0) 
are identical in the finite temperature classical and their 
corresponding quantum dual systems. 



III. GLASSY QUANTUM DYNAMICS 

Response functions in materials of changing plastic- 
ity, e.g., glass formers which as their temperature is 
lowered become progressively more viscous (as well as 
various electronic systems), suggest the presence of a 
distribution of local relaxation times tha t lead to, e.g., 
the canonical Cole-Cole or Cole-DavidsorPSEH and sim- 
ilar forms as we briefly elaborate on. The response of 
a single attenuated mode to an initial impulse at time 
t = scales as g S i ng ie ~ exp(— t/r) with r the relax- 
ation time; the Fourier transform of this response reads 
9singie{<-<j) = ,9o/(l — iw). In systems with a distribution 
(/(r')) of relaxation events, the response functions are 
given by J dr' f(r') exp(— t/r'). Empirically, in dissipa- 
tive plastic systems, relaxations scale as (exp[— (t/r) a ]) 
with a power < a < 1 that leads to a "stretching" of the 
response function. This stretched exponential and other 
similar forms -such as the Cole-Cole (CC) and Davidson- 
Cole (DC) functions- capt ure qu intessence of the distri- 
bution of relaxation times PSES With g(u>) — g G(uj), 
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where go is a constant, the CC and CD forms correspond 
to different choices of G, 

G cc ( w ) = l/[1-M»], 
Gdc{u) = 1/[1 -iur] b . (10) 

The parameters a and b qualitatively emulate the real- 
time stretching exponent c. This distribution of relax- 
ation times might be associated with different local dy- 
namics (dynamical heterogeneities) to which we will turn 
to shortly in subsection |III A| Fluids may fall out of equi- 
librium when cooled at a sufficiently rapid rate (so-called 
"super-cooling") to become "glassy" at sufficiently low 
temperature. As liquids are supercooled, their charac- 
teristic relaxation times and viscosity may increase dra- 
matically. There are several time scales that govern the 
dynamics of supercooled liquids. The so-called "a relax- 
ation" is associated with cooperative motion and leads 
to a pronounced rise of the viscosity (especially so in 
the "fragile" glass- formers). Empirically, in classical su- 
percooled liquids at a temperature T c i, the a relaxation 
times follow the Vogel-Fulcher-Tammann forrrP3 



r(T cl ) 



r e 

DC 



A/(T al -T 



for T > T , 
for T < T . 



(11) 



Here, To is the temperature of the classical system at 
which the relaxation times (if Eq. (Ill is precise) will 
truly diverge and A is an energy scale. Mode coupling 
theory^ and numerous other theories might similarly 
capture aspects of the increase in the a relaxation time. 
By a trivial application of our result of Eq. Q, all of 
these finite temperature classical forms have quantum 
zero-temperature counterparts. Slower dynamics also 
appear as the system density increases. Several works, 
e.gpSl suggest that relaxation times are a function of a 
composite quantity involving both density and tempera- 
ture. It should be noted that all of our derivations start 
from Eq. ([3| . When examining the quantum dual to em- 
pirical forms describing classical liquids, the bare viscos- 
ity of the ambient liquid r\ appearing in these equations 
of motion may, in principle, be allowed to change as the 
temperature (and density) are varied. These may appear 
in addition to changes in T c ; and Vn (capturing, e.g., 
changes in the density). In classical simulated liquids, 
the bare viscosity may be kept constant. 



A. Quantum dynamical heterogeneities 

Our key result of Eq. ([7]) provides a natural bridge 
between classical and quantum correlations between dis- 
parate operators {0(t)}. The remainder of this work will 
largely focus on consequences of this relation. From it, 
correlation times {r a } that appear in the classical prob- 
lem will rear their head in the quantum problem (and 
vice versa). 

We now focus on an intriguing aspect of classical 
glasses which by virtue of the relation of Eq. ^ [as al- 
luded to irf^ leads to the appearance of new dynamical 



correlations in quantum systems. Disorder free models 
for classical glass formers (including various simulated 
quenched systems such as those endowed with various 
classical potentials Vn (e.g.pj) that do not permit simple 
crystalline orders) are known to exhibit "dynamical het- 
erogeneities" (DH)- a non-uniform distribution of local 
velocities® with the location of the more rapidly moving 
particles changing with time. By invoking Eq. ([7|, we 
see that Quantum Dynamical Heterogeneities (QDH) 11 
appear in their corresponding zero temperature quantum 
counterparts. That is, in disorder free quantum systems 
derived (via Eq. ([5])) from the corresponding classical 
systems, zero point dynamics is spatially non-uniform. 

The presence of DH is seen by numerous probes®. One 
often used metric is that of the four-point correlations of 
Eq. j2j in various guises. These correlation functions are 
of the form of Eq. ^ with 0(t) denoting the overlap 
between fields cf> when these arc separated in time, 

o&) = M^-M - (Mt)) <4>-<f(o)}, (12) 



with q any wave-vector. When Eq. ( 12 ) is substituted 
into Eq. ^ , the Fourier space correlation functions (de- 
noted 54 lasslcal (q, t) below) typically have an Ornstein- 
Zernicke type or similar related forms, e.g.,® 



-^classical 



(13) 



with the length scale ^(t) representing the size of the 
typical dynamical heterogeneities when the system is ex- 
amined at two times separated by an interval t. The 
four-point susceptibility is simply Xi(t) = j dx G4,(x,t). 
The key feature of Eq. (131 is that all of the q depen- 



dence has been relegated to a Lorcntzian form while X4 
and £4 are otherwise general time dependent functions. 
We may next invoke Eq. Q to generate the quantum 
counter-part of Eq. |l3]) (or of any other related form) 
and Fourier transform to real space to obtain, in the no- 
tation of Eq. ([2]), the spatial bosonic correlation function 
Q^ u anturn (x—y, t) associated with the potential VQ U antum of 



Eq. ( A12[ ) . The Fourier integral will be dominated by mo- 
mentum space poles at q — ±i£^ . It is clear that in em- 
ploying the transformation of Eq. G 4 5uantum (x - y,t) 
will exhibit exponential decay with the very same correla- 
tion length £4 that is present in the classical system. This 
affords a direct proof of the dynamical length scales £4 
in all zero temperature quantum counter-parts (given by 
Eq. ([5])) to any dissipative classical system that is known 
to exhibit these (and there are numerous known classical 
systems that exhibit dynamical heterogeneities®!). 

Following our mapping, an exponentially decaying real 
space 4-point correlator (r = \x—y\) in the classical prob- 
lem (exp(— r/^4)) will lead to an oscillatory decay in the 
quantum dual. Specifically, if the four-point correlation 
length in the classical problem diverges as £4 ~ t x / z then 
in the corresponding quantum dual system the correla- 
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tion function of Eq. ([2| will be given by 

G 4 ,Quantjm(r, i)~ sin(— r -j- sin — ; 

T X ' Z AZ 



x exp[ 



r tt 

■ —r-r cos — . 
t x I z 2z l 



(14) 



B. Rapidly increasing time scale with concomitant 
slowly increasing length scales in quantum glasses 

There is a proof that a growing static length scale 
must accompany the diverging relaxation times of glass 
transitions^. Some evidence has indeed been found for 
growing correlation lengths (static and those describing 
dynamic inhomogeneitiespSHSS As we noted earlier in 
this work, correlat ion l engths were studied via "point- 
to-set" correlations^^] anc j pattern repetition size^l. 
Other current common methods of characterizing struc- 
tures include (a) Voronoi polyhedrapDHB (b) Honeycutt- 
Andersen indices® , and (c) bond orientation^ all cen- 
tering on an atom or a given link. More recent approaches 
include graph theoretical toola^Sl. Not withstanding cur- 
rent progress, it is fair to say that currently most "natu- 
ral" textbook type length scales do not increase as dra- 
matically as the relaxation time does when a liquid is 
supercooled and becomes a glass. 

It is worthwhile to highlight that one of the most per- 
tinent naturally increasing length scales is that associ- 
ated with the typical size of the dynamical heterogene- 
ity (i.e., £4 of Eq. (13)). Similar to other measures, 
this typical length scale does not increase as rapidly 
as the characteristic relaxation time does as the glass 
transition is approached. Recent work for a three- 
dimensional bi-disperse repulsive glass^ [with the pair 
potential V^(|a?|) = e(a a b/r) 12 between two particles (a, b) 
of two possible types ((a, b) G 1, 2) with cr ab = (a a + cr b )/2 
and oil<J\ — 1.2)], suggests that 



exp(fc£|), 



(15) 



with ~ 1.3 and k a constant. (An alternate assumed 
algebraic form r <~ £| leads to an extremely large dy- 
namical exponent z ~ 10.8.) In these cases, the dynam- 
ics changes dramatically with little notable change in the 
spatial correlation length. 



C. Quantum Critical Jamming 

The mapping between dissipative classical and quan- 
tum systems also suggests a new quantum critical point 
in related systems. The classical jamming transitiorPSHSS 
of hard spheres/disks from a jammed system at high den- 
sity to an unjammed one with spatially heterogeneous 
motion at lower densities is a continuous transition with 
known critical exponents, both staticP^'and dynamic^. 
In most solids, the ratio of the shear to bulk modulus is a 
large number like 1/3. However at the jamming thresh- 
old, this ratio tends to zero. This renders the system very 



soft to shear stresses. This is one of the things that makes 
jamming different from other solids^. By our mapping, 
the classical jamming transition has a quantum analog. 
Replicating the mapping of the previous subsection (and, 
in particular, Eq. |7]) therein), wc may derive an analog 
quantum system harboring a zero temperature transition 
with similar critical exponents. The c lassical zero tem- 
perature critical point ( "point J" )) 48 l 5 °l may rear its head 
anew in the form of quantum critical jamming (at a new 
critical point — "Quantum point J") in bosonic systems. 
A schematic of the phase diagram of the associated quan- 
tum system is depicted in Fig. [T] 

We may ascertain dynamical exponents from those re- 
ported for the classical jamming systeirP^. The classical 
low temperature temperature system (T c i 0) maps, 
according to Eq. ([5]), onto a zero temperature quantum 
system in its large mass limit. Bosons of infinite mass are 
not trivial due their statistics. Specifically, for a classi- 
cal system of mono-disperse soft spheres with a repulsive 
force that is linear in the amount of compression, it was 
found that the correlation length £ and relaxation time 
r scale 59 as 



(16) 



In Eq. ( 16 ), p denotes the density with pj being the criti- 
cal density at the jamming transition. Eq. ( 16 1 describes 



how the spatial and time scales diverge as the density is 
increased and approaches (from below) the density at 
the jamming transition. The correlation length in the 
jammed systems is set by the scale at which the number 
of surface zero modes is balanced by bulk effects. Taken 
together, these imply that, on approaching the transition, 
the relaxation time increases much more rapidly than the 
correlation length, r ~ £f with a very large effective dy- 
namical exponent z ~ 4.6. By use of Eq. 0, the same 
behavior is to be expected for the quantum system gov- 
erned by the corresponding quantum potential VQuantum- 
The classical (and thus quantum) jamming exponents are 
the same in two and three-dimensions. It may be re- 
marked that a similar dynamical exponent was found for 
a Bose glass model suggested to describe vortex line in 
high temperature superconductors^. In physical terms, 
for charged bosons, the jamming transition constitutes 
a transition from a metallic system (when the system is 
unjammed and behaves as a fluid) to a jammed state (an 
insulator). We note earlier work rationalizing metal to 
insulator transitions in terms of electron pinning^. In 
the bosonic jamming that we describe above, no pinning 
is present and the transition is driven by particle inter- 
actions. 

As in earlier sections, we see that time scales increase 
far more precipitously than spatial correlation lengths. 

One of the hallmarks of jammed systems is that the 
spectral density of vibrational excita tions, D(u>) is con- 
stant at the jamming threshold 5 ^. (I n a Debye solid 



D(u) 



ijj- 



.) This is independent of potential, dimension 



of the system and size of the system. As one goes away 



() 



T H ~ 1/mass 




1 /density 



FIG. 1. (Color online.) The phase diagram of the zero tem- 
perature quantum jamming transition with line of J points. 
The phase boundaries and axis were formed by employing 
the phase diagram of the classical syste and examining 
the duality between the classical and quantum system- i.e., 
comparing the parameters in the quantum sys tem of Eq. |5| 
with the classical system defined by Eqs. ( 3|4 1. 



from the jamming threshold, D(uj) has a plateau down to 
a frequency lo* . Below u* , the density drops and starts to 
look more like a Debye solid. This frequency u>* veers to 
zero on approaching the transition. The number of low 
energy modes is set by the absence of constraints and 
Maxwell counting arguments. This constant density of 
states implies an enormous increase in the low frequency 
excitations. These can be measured by looking at the 
trajectory of a single particle and Fourier transforming 
its motion. One can also look at the spatial extent of 
these modes. They too are different from what one nor- 
mally gets. They are quasi-localized (or resonant) at low 
frequency below u* . Above u* they are extended but 
still do not look at all like plane wavesP^ The results 
concerning the density of the modes hold unchanged for 
the dual quantum system. This is so as the dynamical 
matrix D is formed by the second derivatives of the po- 
tential Vn relative to the displacements. Thus, the same 



statement about mode density of states that appear in 
classical jammed system is reflected verbatim in the clas- 
sical potential Vpjf which we use to construct the quantum 
many body potential VQ ua ntum({^}) from Eq. (J5j|. Any 
appearance of zero energy (bulk or surface) modes in the 
classical system will identically hold also in the quantum 
system. 



IV. LATTICE SYSTEMS 

Thus far, we focused on continuum viscous classical 
systems which, as we have seen, mapped onto continuum 
bosonic systems. We briefly remark here on classical lat- 
tice systems which similarly exhibit dynamical hetero- 
geneities and a jamming type transition. [In appendix 
|H| we further brief ly exp and on our derivations for lat- 
tice systems. ] RefsPMl studied, respectively, the 2DN3 
and 3DN2 models on the square and cubic lattice models 
in d = 2 and d = 3 dimensions. In the square lattice N3 
model, particles are endowed with hard core repulsive in- 
teractions that extend up to a distance of three steps on 
the lattice. Similarly, in the cubic lattice N2 model parti- 
cles cannot be nearest neighbors nor next nearest neigh- 
bors (i.e., the repulsive hard-core interactions extend up 
to a distance of two steps on the lattice). Thus, simi- 
lar to the continuum systems that we discussed earlier, 
these models may be regarded as that of classical hard 
core spheres. Following the mapping reviewed in section 
|lH the quantum counterpart of such systems is that of 
dominant hard sphere interactions augmented by con- 
tact stick y inte ractions. In the classical systems, simula- 
tion start ^£4165| w it n an infinitely fast quenching wherein 
particles are added whenever possible and diffuse other- 
wise; this process is halted when the desired density is 
reached. A clear increase was noted in the le ngth scales 
that characterize the dynamical heterogeneity^ 6 ^. The 
continuum jamming transition discussed earlier may have 
a lattice counterpart for Cooper pairs as we now elabo- 
rate on. A natural quantum counterpart to th e N3 (N2) 
model is given by an extended Bosc Hubbar d 66 * 67 * type 
model with infinite hard core repulsions, 



H 



-t^ibjbj +h.c.) + U^n l {n l 



1) 



(ij) 



(17) 



where Vij — > oo for lattice sites i and j which are fewer 
than four (or three) steps apart and the onsite Hubbard 
repulsion U is divergent (U — > oo) as well. The Hubbard 
term leads to a penalty only when there is a double or 
higher occupancy. Based on our considerations thus far, 
we expect to obtain the quantum bosonic counterpart to 
the classical jamming transitions found in the classical 
2DN3 and 3DN2 models. This bosonic system may have 
all of the characteristics of the classical jammed system 
including dynamical heterogeneities and a large dynami- 
cal exponent z. For completeness, we briefly comment on 
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the difference between the lattice system of Eq. (17) and 
the "Bose glass" first introduced irPH The Bose glass ap- 
pears in the bare (i.e., that with = 0) disordered ren- 
dition of Eq. (171 with the general Bose Hubbard Hamil- 



tonian (with general finite repulsion U) being further 
augmented by a local chemical potential term — \iirii 
wherein fa is a spatially non-uniform random quantity. 
By contrast, the lattice Hamiltonian of Eq. (17) as well 



as the continuum models that we discussed in earlier sec- 
tions are free of disorder. The amorphous characteristics 
that these clean systems may exhibit arc borne out of 
"self-generated" randomness^ — not randomness that is 
present in the parameters defining the system. 



V. ELECTRONIC SYSTEMS WITH PAIRING 
INTERACTIONS 

Up to now, building on and extending the mapping 
between classical dissipative systems and zero tempera- 
ture bosonic theories, we focused on hard core bosons. 
We now turn to the ground states of Fermi systems. In 
particular, in this section, we will consider standard elec- 
tronic systems with pairing interactions, 



H 



ka C ka 



V k,i C kx c - k 'i c -Tl c Tv 

k,l 



(18) 



where a =t,| is the spin polarization index and Vki is 
the interaction strength between the Cooper pairs 
\k f;—k l) and \l X) ■ As is well known (and is 

readily verified), the following Fermi billinears 



* Art -fej. 



(19) 



corresponding to the creation/annihilation of Cooper 
pairs satisfy hard core Bose algebra. We next consider 
what occurs if, within the ground state, the occupancies 
of the single particle states are correlated inasmuch as the 
electronic states on which the standard pairing Hamilto- 
nian of Eq. ( 19 ) operates can be created by applications 



of Cooper pair creation operators on the vacuum (i.e., if 
the ground state is invariant under the combined opera- 
tions of parity (k — > —k) and time reversal (a — > —a))- 
When the ground state is strictly invariant under the 
combined effect of these symmetries, we may express the 
Hamiltonian of Eq. ( 18 ) as a bilinear in the hard core 
Bose operators, 



H 



(20) 



The hard core (Fourier space) Bose algebra of the cre- 
ation and annihilation operators [as, in particular, man- 
ifest in the relation (&g) 2 = mandating that no more 
than one boson can occupy any given (Fourier space) 



site] is identical to that of raising and lowering opera- 
tors in the spin S = 1/2 system. Thus, a simple ex- 
tension of the standard real-space Matsubara-Matsuda 
transformational is given by 
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(21) 



Substituting Eq. |21J) into Eq. (20), we arrive at an 
XY model. In situations in which the band dispersion 
is nearly flat (and may be omitted for fixed particle 
number), in determining the ground state(s), we must 
only find the pairing V that affects pair hopping. Similar 
considerations apply in real-space when Cooper pairs are 
short ranged and may be replaced by real-space hard- 
core bosons. Hard core real space contact interactions 
correspond to uniform Vt f (independent of k and I) as in 
the BCS form for the pairing interactions. In such cases, 
whenever the system is dominated by hard core contact 
interactions between the bosonic Cooper pairs we see, 
replicating our analysis thus far, at zero temperature, 
that the system may undergo a jamming type transition 
between an itinerant and jammed phase at sufficiently 
high densities or pressure. In this case, it displays rapidly 
increasing relaxation times concomitant with spatial cor- 
relations that do not increase as dramatically as the re- 
laxation times do on approaching this transition. 



VI. GENERAL BOSONIC QUANTUM 
CRITICAL POINTS WITH CLASSICAL 
ANALOGS 

Up to this point our focus has been on supercooled 
systems. Eq. Q implies that also standard (non glassy 
or spin-glassy) classical transitions have corresponding 
zero temperature quantum analogs. Bona fide critical 
or discontinuous classical transitions have quantum du- 
als according to the methods that we invoke. Ising type 
transitions have critical exponents z which are given by 
Z2D — 2 and a somewhat smaller value in D = 3 dimen- 
sions. At the critical point of the liquid to gas transition, 
the system is in the Ising universality class. Thus, the 
transition from a quantum Bose gas to a fluid is governed 
by the same dynamical exponents. 



VII. POSSIBLE IMPLICATIONS FOR 
EXPERIMENTAL DATA 

We now, very briefly, turn to a discussion of possible 
data analysis of experiments. One of the main messages 
of our work is that classical physics associated with over- 
damped classical systems can rear its head in the quan- 
tum arena. Correspondingly, data analysis which has 
led to much insight in the study of classical glasses and 
other damped systems may be performed anew for quan- 
tum systems. A principal correlation function which we 



focused on in this work has been that of the four-point 
correlation function of Eq. ([2| . This correlation function 
need not be directly measured in real-time. For instance, 
scanning tunneling spectroscopy (STS) data taken at dif- 
ferent positions and bias voltages may provide a valuable 
conduit towards the evaluation of the four-point correla- 
tor when it is expressed as an integral over frequencies (or 
associated bias voltages). Rather trivially with </>(x, V) 
denoting the local density of states at location x for a 
bias voltage V, and e* the electronic charge, the corre- 
sponding four-point correlation function is given by 

G 4 (x-y,t) = J dV 1 dV 2 dV z dV^ 

(5ct>(x, V 1 )dc/>(x, V 2 )Sm V 3 )Sc/>(y, T/ 4 ))e (e * t(yi+V3) 

-C 2 (<),(22) 

with the two-point auto-correlation function 

C(t)= [ d d x [ dV(Scj)(x,V)S(j)(x,0))e iettv . (23) 



As we have discussed in Section |III A| quantum dynam- 
ical heterogeneities may be manifest in this correlation 
function via the standard relation of Eq. (13). 



VIII. CONCLUSIONS 

A central result of this work is the exact temporal cor- 
respondence of Eq. ([7]) that spans both equilibrium and 
non-equilibrium dynamics in general time dependent sys- 
tems (so long as either the initial or the final state of the 
classical system is that of thermal equilibrium). This 
equality relates (i) the auto-correlation function of Eq. 
([6]), for any quantity O when evaluated for the classical 
dissipative system of Eq. Q with a many body poten- 
tial energy Vjv, to (ii) the auto-correlation function of 
the very same corresponding quantum operator O in a 
dual bosonic system governed by the Hamiltonian of Eq. 
([5]) . When fused with known results for dissipative classi- 
cal systems, this extremely general equality immediately 
leads to numerous non-trivial effects which we introduced 
and readily proved as a matter of principle. These in- 
clude: 

• Quantum dynamical heterogeneities (QDH). We il- 
lustrated that similar to classical systems even in the 
absence of disorder, bosonic systems can, at zero tem- 
perature, exhibit spatially non-uniform zero-point mo- 
tion. Of course, in translationally invariant systems, the 
average (time averaged) dynamics is uniform. However, 
at any given time, there are particles that move more 
rapidly than others. We suggested how experimental 
data may be analyzed to search for quantum dynamical 
heterogeneities in electronic systems. 

• The length scale characterizing the zero temperature 
QDH, the four-point correlation length £4 (a trivial ana- 
log of its classical counterpart) may increase as the dy- 
namics of the clean Bose system becomes progressively 



sluggish. However, albeit its rise, this length scale may 
increase much more slowly than the relaxation time. The 
far more rapid increase of the relaxation time as com- 
pared to readily measured length scales is a hallmark 
of many electronic systems. Cast in terms of quantum 
critical scaling (if and when it might be realized) , the ef- 
fective dynamical exponent z capturing the relation be- 
tween correlation lengths and times is very large [z 3> 1). 
Other relations such as those of Eq. ( 15 ) may hold once 



they are established for viscous classical systems. 

• Similar to classical systems, quantum systems may 
jam at high densities or pressure notwithstanding zero 
point motion. The character of the jamming transition 
in zero temperature quantum systems is identical to that 
of that of their corresponding classical finite temperature 
counterparts. As the classical systems exhibit a critical 
point at the jamming transition (at "point J") so do their 
bosonic counterparts. As a result, we established the ex- 
istence of a new quantum critical point — associated with 
a quantum critical jamming of a hard core Bose system. 
As in the other systems that we discussed, the character- 
istic relaxation time diverges more precipitously than the 
correlation length on approaching the transition ( "Quan- 
tum point J" ) with a large effective dynamical exponent 
2-4.6. 

• The continuum theories that we predominantly fo- 
cused on may have a broad applicability as continuum 
theories describe the same physics as their lattice rendi- 



tions do in the vicinity of critical points. In Section IV 
we discussed specific possible lattice renditions. 

• The results that we derived for zero temperature 
bosonic theories suggest similar features in electronic sys- 
tems. In some cases, as discussed in Section |Vj finding 
the ground states of interacting electronic systems can be 
cast in terms of a corresponding zero temperature hard- 
core Bose problem. 



• The general mapping of Eq. (21) between hard core 
bosons to S = 1/2 spin systems (in either momentum (k) 
or real (r) space) along with the complementary relation 
for the z component of the spin, 



fish 



1/2] SI 



(24) 



in the same space, allows us to derive similar results for 
certain spin S = 1/2 systems. 

Thus, with the aid of the viscous classical-many body 
quantum correspondence of Eq. Q, we trivially estab- 
lished all of these results without the need to perform 
various standard and far more laborious computations 
for quantum systems. 

Other possible extensions of our results include quan- 
tum turbulence from known results on the dynamics in 
classical fluids that exhibit turbulence, the relation be- 
tween localization (or caging) in classical systems and 
their corresponding quantum counterparts. We may sim- 
ilarly examine disordered systems- for a random classi- 
cal potential Vn, the corresponding quantum potential 
VQuantum is also random. 
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As we will further detail in the appendix (see Eq. ing position vectors {x(t)} at time t is given by 



(A20)), if we are given a known quantum ground state, 



we may find the corresponding effective classical poten- 
tial. With the aid of calculations on how the correlation 
functions of the classical system depend on time as pa- 
rameters in the classical potential are varied, we may 
then determine the corresponding time dependent corre- 
lation functions of the many quantum system. That is, 
we need not always find corresponding quantum systems 
to classical systems; the Fokker-Planck mapping also en- 
ables to go in the opposite direction from quantum sys- 
tems to classical ones. 

Numerous related extensions may be considered. For 
instance, we may consider magnetic and other systems in 
which fcrmionic degrees can be formally integrated out 
leaving only effective bosonic degrees of freedom. Con- 
sequences for the Ward identity relating four-point with 
two-point correlation functions may be considered. The 
Langevin equation may be re-examined for single vortex 
crossing of a narrow superconducting wire at finite tem- 
perature to derive the mapping for the quantum dual at 
absolute zero temperature! 73 * 7 ^ This would offer an alter- 
native path for exploring the viability of quantum phase 
slips in nanowires. 
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Appendix A: The Fokker Plank of the dissipative 
classical system and its relation to a quantum 
Hamiltonian 



In w hat fo llows, we highlight features of the 
mapping * 15 ! 16 * 20 ! between classical dissipative systems 
and bosonic systems. We first set the preliminaries 
following 16 . Given the initial vector xo of the coordi- 
nates of all particles at time t = f , the time dependent 
probability distribution V(x, t; Xo,to) for the correspond- 



N 



7>({f},i;{f },i ) = fll *[{*»•(*) -*i}]>W,{2o}J(Al) 

i=l 

where (— )n,{x } denotes the average over the random 
noise r\ (which we will take to be the Gaussian white 
noise of Eq. Q) given that initially, at time t = t Q , the 
particle coordinates were {xq}. The average of a general 
function 0({x(t)}) is then 



d dN xV({x},t;{x },t )Om) 
= (0({f(t)})) Wl{Jo} . 



(A2) 



It is convenient to write V({x}, t; {afo}, to) in a Dirac 
notation as 

V({x},t;{x },t ) = ({x}\P(t,t )\{x }). (A3) 

Time translation invariance and the Markov property of 
these probabilities, 

d diV : r'({ : r}|P(t,0|{5'})({^}^(t , J to)|{^o}) 

= / d m x'V({x},t;{x'},t')Va^},t';{^},t ) 

= V{{x},t;{x },t ),(A4) 

imply that 



P(t,t )=Te 



- P t H FP {t')dt' 



(A5) 



with a "Fokker-Planck" Hamiltonian Hpp and where T 
is the time ordering operator. We now return to our 
particular classical to quantum mapping. The summary 
below closely follows this mapping as presented by Biroli 
et alP^l. In what follows, we set 



P = P(t,to)\{x }). 



(A6) 



For the classical dissipative system of Eq. (J3|) , the 
probability distribution P({x}) evolves according to the 
Fokker-Planck equation 



dP 
~dt 



H FP P, 



where the Fokker Planck operator is 



Hpp 



1 



V 7i 



ViVN +T cJ V, 



(A7) 



(A8) 



with T c i being the temperature of the classical system 
[setting the noise strength in Eq. Q]. Eq. (A8) follows 
from a direct differentiation of Eq. ([AT]) while invok- 
ing Eq. ^ for the derivatives of the coordinates {xi(t)} 
in the argument of the delta functions and performing 
short time averages. (Thus, this equation and our re- 
sults pertain to systems in which the dynamics is suffi- 
ciently slow such that short time averages over the noise 
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•q at fixed temperature are sensible.) A detailed deriva- 
tion of the Fokker-Planck equation for this and more 
general Langevin processes appears in many excellent 
textbooks, e.g.j^l. The operator Hpp is non-Hermitian. 
Each eigenvalue is generally associated with differing left 
and right eigenvectors. The Fokker-Planck equation can 
be mapped into a Hermitian Hamiltonian byS21 



H = e^'^H 



F pe 



-Vn/(2T c! ) 



(A9) 



if the second derivatives are exchangeable, V^Vj-Vn = 
VjViVNpD A direct substitution leads to the quantum 
many body Hamiltonian of Eq. Note that, thus 

far, we have allowed Vn to be completely general. This 
potential energy may include one body interactions (i.e., 
coupling to an external source) , pair interactions between 
particles, and three- and higher-order particle interac- 
tions. A key point that we will further invoke later is that 
the transformation of Eq. ( A9 ) leading to a Hermitian 



quantum Hamiltonian can be trivially performed at any 
given time slice when Vfj and T c i are, generally, time de- 
pendent. It is worth highlighting that in non-equilibrium 
time dependent classical systems, the temperature T c i is 
set by the time dependent noise amplitude (following Eq. 
Q). We consider what specifically occurs when the clas- 
sical potential energy in Eq. ([3| is the sum of pairwise 
interactions, 



Vn ({*}) = \zZ V ^ 3 



Xj). 



(A10) 



For such systems, the quantum many body Hamiltonian 
of Eq. ([5]) explicitly contains an effective potential which 
is the sum of two and three body interactions, 



uantum <{*}) 



E 1 



N 



4T r 



-(V^ N ) 5 



-2g^ + S k«, ■ <A11) 

For a given classical two body potential in d dimensions 
which is both translationally and rotationally invariant, 
V(x) = V(\x*\), the resulting quantum potential energy 
is given by (as irP^l yet now trivially extended to include 
general T cl ), 



VQ U antum({a ? }) — „ £ ^Quantum ( 



, \ ^ 3 — body / -> -* -> -> \ 
' 2-^1 ^Quantuml^i ~ X h X 'i ~ X j' h 

^(f)^V 2 ni) + ^r[W(f)f 

= - — V'{r)-V"{r) + ^-[V'{r)f- 
r 2T ci 



v 3 - body (x,x') = — V7(a;) • VV(x' ) 



1 



4T ci 
1 x x' 



7 V'(r)V'(r>), 



(A12) 



with r = |a?| . For a classical potential V(r) = 
Vq exp(— \[{r/a) 2 — 1]) the corresponding pair term in 
Eq. (|A12|) is 



'" i "" r (r) = 



-V(r) 



T r] (J 



In the limit A — > 00, the classical system corresponds to 
that of hard spheres where a is the diameter of the hard 
sphere and the quantum potential of Eq. ( A13 ) similarly 



exhibits a dominant hard sphere repulsion (augmented 
by an attractive potential at the sphere boundaries that 
is of range 1/A). 

For gen eral Vn({x}), the Fokker-Pla nck o perator of 
Eq. (A8) has non- negative eigenvalue d 15 * 16 l For any 
time independent Hpp, the zero eigenvalue state — i.e., 
the ground state — which according to Eq. (jATb corre- 



sponds to a stationary (time independent) probability 
distribution P. This is the equilibrium Boltzmann dis- 
tribution 



-I3V H {{3}) 



(A14) 



with Zjq the partition function associated with Vn ({2}) 
and P = 1 /T c i . This is readily rationalized by the follow- 
ing argument. For a finite size system, the linear eigen- 
value equation 



FPJbc 



(A15) 



with the matrix row/column indices b and c denoting 
classical configurations, has a (finite size) matrix Hpp 
with positive off-diagonal elements and negative diagonal 
entries. Specifically in Eqs. ( A7|A15| ), the transition 
matrix Hpp has entries that relate the probabilities of 
going from state b to state c in a given (infinitesimal) 
time interval. If these states are different (b 7^ c) then 
clearly (Hpp)b c > 0. The diagonal elements (Hpp)bb 
provide the probabilities of "leaking out" of state b and 
going to all other states c^b. From all of this it follows 
that 



{Hpp)bb 



b'^b 



FPjbb 



)bb> < 0. 



(A16) 



Detailed balance, i.e., the fact that the probability of 
going from b to c is the same as that of going from c to 
b, asserts that 



(Hpp)bce 



(Hpp) c i,e 



-/3-Ec 



(A17) 



Eqs. (A15 A17) illustrate that the Hamiltonian of Eq. 
(A9) is Hermitian. In this classical system of Eq. pj), 
the energies of the classical states E c are simply given 
by Vn({x}) evaluate d for the c lassical configurations c. 
With the aid of Eqs. ( A16 A17 ), it is easy to see that the 



column vector p^i ml = Z N exp(—f3E c ) (i.e., the distri- 



bution of Eq. (A14)) is a null eigenvector of Eq. ( A15 1 



This probability eigenvector corresponds, of course, to 
the equilibrium Boltzmann distribution. The factor of 
ZZ- is inserted to ensure normalization of the classical 
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probabilities (for any eigenstate): ^2 C P C = 1- Now, we 
can add a constant to the finite dimensional matrix 



Hpp — > Hpp — const. = H 



(A18) 



to generate a matrix {—H FP ) that has all of its elements 
positi ve (—H FP )bc > 0. Specifically, to this end, in Eq. 
(A18 we can choose const, to be any constant larger 



than the sign inverted smallest off-diagonal element of 
(—Hpp), i.e., const. > — mhi^^Hpp}^. For such a 
positive matrix, we can apply the Perron-Frobenius the- 
orem which states that the largest eigenvector of (—H FP ) 
is non-degenerate and that eigenvector is the only eigen- 
vector that has all of its elements positive with all other 
orthogonal eigenvectors having at least one negative ele- 
ment. Clearly, all of the eigenvectors of Hpp and H FP 
are identical with the corresponding eigenvectors of both 
operators merely shifted uniformly by a constant. With 
all of the above in tow, we see that p ec i ml corresponds 
to the largest eigenvector of (—H FP ) and is thus also the 
largest eigenvector of (—Hpp). For a time independent 
Hpp, as p e i ul1 was the null eigenvector of Hpp, it fol- 



lows that all other eigenvalues of Eq. (A15), e > 0, are 
positive and, according to Eq. (A7) and explicit earlier 
discussions evolve with time as exp(— et) —^t-^oo 0. Thus, 
physically (as it to be expected) at long times the sys- 
tem attains its equilibrium configuration of p e i utl . In 
the corresponding zero-temperature quantum problem, 
the dominant classical equilibrium state with a lowest 
energy. We will thus label it in Appendix [Cj by \G). The 
transformation of Eq. ( A9 ) relates the operators in the 



classical Fokker-Planck and zero temperature quantum 
problem to one another. The transformation for the right 
eigenvectors P of Hpp, which we explicitly denote below 
as \—)fPi t° the eigenvectors of the quantum Hamilto- 
nian H is trivially 



)FP 



exp(-lV(2T c( )))|- 



)FP 



Similarly, the left eigenvectors ((— \fp) of Hpp are to be 
multiplied by exp(Vi\i/(2T c j))) in order to pass to the left 
eigenstates of the quantum problem. (In the quantum 
problem defined by H, the left and right eigenstates are 
trivially related to each other by Hermitian conjugation.) 
Applying Eq. (A19) to the null right eigenstate of the 



Fokker Planck Hamiltonian of Eq. ( A14 1, we see that the 



quantum eigenstate of H corresponding to this classical 
equilibrium state is given by 



*o(M) = ^^cxp(- 

V 



1 

2T r 



V N ({x})). (A20) 



The prefactor in Eq. (A20) is set by the normalization 
of this quantum state. When comparing Eqs. ( A14|A20 ) 
to one another, we see that this wavefunction is related 
to the classical equilibrium probability eigenvector by 
the appealing relation ^({cc}) = y/ P e Q uil ({x}) When 
Vn ({:?})) is symmetric under the interchange of parti- 
cle coordinates, the resulting wavefunction may describe 
bosons. For two body interactions Vij (Eq. AlOl that 



are symmetric under the permutations of i with j, the 
wavefunction of Eq. ( A20 1 is symmetric under any per- 



mutations of th e part icles. Thus, the ground state wave- 
function of Eq. ( A20 ) is a Jastrow type wavefunction de- 



scribing a bosonic system. Of course, generally, Vn({x})) 
can include not only two body terms but also single body 
contributions (local chemical potentials or fields) as well 
as three- and higher-body interactions. Although ob- 
vious, it is worth noting that if Vn ({:?})) (and thus the 
quantum Hamiltonian H) is invariant under any pairwise 
permutation Py, i.e., if [H, Pij] = then the symmetry 
of the initial wavefunction \Vq (corresponding to the clas- 
sical Boltzmann distribution for a system initially at an 
equilibrium at temperature T c i) does not change as the 
system evolves with time (including general arbitrary H 
corresponding to classical variations in temperature and 
other parameters). 

We complete this appendix with a brief discussion con- 
cerning a generalization of Eq. ( A9 1 . It is, of course, pos- 



sible to write down a general similarity transformation, 



H 1 — S ^HppS 



(A21) 



with a time dependent operator S. Under Eq. (A21|, 
the Fokker-Planck equation of Eq. ( A7 ) reads 



ft* 



-if'*, 



(A22) 



where \T> = S 1 P with P given by Eq. (A6) 



Appendix B: Supersymmetry and the mapping from 
classical stochastic systems to quantum mechanics 

As is well known, there exists a beautiful link between 
stochastic classical statistical mechanics and supersym- 
metric quantum systems, e.g.p- 2 . This connection is es- 
pecially immediate for the ground state wavefu nctio ns 
which are of zero energy (as indeed that of Eq. ( A20[ )). 
This might lead to the impression that the results de- 
rived using the correspondence between classical dissipa- 
tive systems and quantum duals are rather limited and 
special and was informally suggested by some to lead 
to un-normalizable wavefunctions if non-constant equi- 
librium classical states are considered. As we will ex- 
plicitly elaborate on in Appendix [Cj the time evolution 
operator U(t) in the corresponding dual quantum prob- 
lem is unitary and thus if an initial state is normalized 
(such as that of Eq. ( A20 ) corresponding to an initial 
classical equilibrium state) then the quantum state will 
remain normalized at all positive times [and vice versa 
for a unitary evolution towards a final state of the form 
of Eq. (A20)]. Although the ground state of the quan- 



tum problem may be dismissed as trivial and special, the 
relations concerning the time evolution to states that are 
not of the form of Eq. ( A20 1 are not as immediate. These 



relations concerning the dynamics form the core of this 
work. For pedagogical purposes, we very briefly review 
here some central notions concerning the mapping of the 
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Fokker Planck process to supersymmetric quantum me- 
chanics as they, in particular, pertain to the equilibrium 
problem. The explicit use of supersymmetry will not be 
invoked in the below and the discussion will be made as 
simple as possible. The Hamiltonian of Eq. ^ can, for 
fixed 7, = 7 in a simple (single-particle) one-dimensional 
rendition which we adopt for ease of notation, be written 
as 



7 



(Bl) 



where, 



At 



.4 



d_ V 
dx~ + 27V 

_d_ V 



(B2) 



In the higher-dimensional many body problem, the gra- 
dients are relative to each of the Cartesian coordinates of 
all of the particles and V is replaced by Vn ({&}). Clearly 
jA'A > and thus if a zero-energy eigenstate of H can 
be found it is the ground state. Now, the square root 
of the classical equilibrium distribution function, i.e., the 
wavefunction of Eq. ( A20 1 is clearly a null eigenstate of 
the operator A above. Inserting Eq. (B2| into Eq. (Bl) 



leads to the identification of the quantum many body 
potential in terms of the corresponding classical poten- 
tial energy Vn({x}). The astute reader will, up to triv- 
ial alterations, recognize these operators as the standard 
raising and lowering operators of the harmonic problem 
when V is a harmonic potential. We briefly return to this 



point in subsection F 2 



Appendix C: Derivation of the quantum to classical 
correspondence for general two operator correlators 

The central role of this appendix is the derivation of 
Eqs.( C18|7 | (or, more precisely, the derivation of Eqs. 



(C16 



C17| ) that lead to Eqs. (|7} |C18[ ). The s ole assump- 



tion made in the below derivation of Eqs. (C16 C17 
C18 1 is that the classical system starts from its equi- 



librium state and then evolves with some general (time 
dependent) potential Vn(£). This will be mapped onto 
analytic continuations of the correlation and response 
functions of a quantum system that starts at time t = 
in its ground state of Eq. ( A20 ) and then evolved with 



the corresponding (time dependent) Hamiltonian H(t). 
It is important to emphasize that we make no assump- 
tions regarding the final (and intermediate) states. The 
classical (quantum) system need not stay in equilibrium 
(or ground state) as it evolves in time. Before detailing 
the derivation, we collect the basic relations discussed in 



Appendix |A"| with several new definitions, 

P({x},t) = {{x}\P(t)), 

e -V N ({x})/T cl 

P G ({x},t) = ({x}\G) = , 

an 

H FP \G) = 0, 

(+IW> = i, 

H = e v »/ 2T °>H FP e- v «/ 2T ". 



(CI) 
(C2) 

(C3) 
(C4) 
(C5) 



These will serve as a point of dep arture for the calcu- 
lations in this appendix. Eq. ( |C1[ ) represents a general 
probability distribution in bra-ket notation. Eq. ( C2 ) de- 



fines the ground state distribution as a Boltzmann distri- 
bution in bra-ket notation. Eq. ( |C3| ) defines the ground 
state as the eigenvector of the Fokker-Planck Hamilto- 
nian with zero eigenvalue. Eq. (C4) defines the state 
|+) to be the uniform state such that |+) = j d{x}\+). 



Lastly, Eq. (C5) can be used to find a relationship be- 
tween H F p and H FP . 

Armed with these, we now proceed to some simple cal- 
culations. As H is Hcrmitian, 



-V N /2T c! u-t „V N /2T C , 



/ft: 

(= H)=e v -/ 2T ^H FP e- v ^'. 



(C6) 



Explicitly multiplying by e v "/ 2Te ' on the left and by 
e -V N /2T cl Qn tne r jg n t i eac i s to 



Ht B = e v »/ T <*H. 



l FP 



FP e 



-v N /T ci 



(C7) 



We will now prove that the state |+) is a left eigenstate 
of the Fokker-Planck Hamiltonian with zero eigenvalue. 
Beginning with a simple extension of the definition of the 
ground state, 



H FP \G) = 0^ (G|/ft p = o. 
As is evident from Eq. |C7|, this is equivalent to 



( G \ e V»/T clHFpe -Vu/T cl 



0, 



which (from Eqs.(C2 C4|) implies that 



J N 



\H FP e 



-Vn/TcI 



0. 



(C8) 



(C9) 



(CIO) 



This illustrates that this uniform state is a left null eigen- 
state for all Fokker-Planck Hamiltonians (i.e., (+\H F p — 
0). 

We will now derive our central result of Eq. Q . To- 
wards this end, we write anew the classical correlation 
function of Eq. ([6]), 



G, 



classical 



(t) = (Oi(t)O 2 (0)>. 



(Gil) 



By Bayes' theorem, the joint probability distribution, 
P({x},{y}) = P({x}\{y})P({y}), the probability of 
finding coordinates {x} at time t and coordinates {y} 
at time 0, is given by the product of the conditional 
probability of finding {x} at time t given {y} at time 
with the probability of attaining {y} at time t = 0. 
For a lattice system with fields <f> at different lattice 
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sites (which we will briefly return to in Appendix 
gj, the equality P({x},{y}) = ^(MIM)P(M) 
is to be replaced by P({(p(x, t)}, {4>(x, 0)}) = 
P{{(j)(x,t)}\{(j)(x,0)})P{{(l)(x,0)}). As discussed in 
Appendix \A» , the ground state has a probability distri- 



-/3Vn(M) 

see Eq. (C2)]. The conditional probability -P({^}|{y}) 



bution given by a Boltzmann distribution Z 



N 



can be expressed in terms of the matrix element of 
Te~ h Hfp ^ ) dt (where T is the time ordering operator) 
as this conditional P satisfies Eq. (A7|. With 0\ 



depending on the coordinates {y} at time t and 2 on 
the coordinates {x} at time t — 0, all this implies the 
form of the expectation value of Eq. ( Cll ), 



G c 



lassical 



(t) 



d^xd^yO x P{{x},{y})0 2 



= d dN x d dN y O x P{{x}\{$})0 2 Pm) 



d dN x d dN y O x {{x}\Te-ti HFp W dt '\{$})0 2 

e -PV N ({y}) 



It is important to re- iterate and emphasize yet again that, 
in the last line above, we only assume that the initial state 
{{y} at time t = 0) is in thermal equilibrium. The system 
need not be in thermal equilibrium at positive times. As 
stated earlier, this is the sole assumption made in this 
derivation for general time dependent systems with dy- 
namical Vn (and thus a time dependent Fokker-Planck 
operators) . A similar derivation would hold mutatis mu- 
tandis when the system is initially out of equilibrium and 
is in eq uilib rium in its final state. 

Eq. (Clk asserts that / d dN x({x}\ = (+|. Invoking 
this along with Eq. (|C2|), we have 



G cla5sical (t) = (+\G 1 Te-ti H ^< it '0 2 \G). (C12) 



As is evident from Eq. (CIO), inserting an exponentia- 
tion of Hpp to the right of the state (+| leads to a mul- 
tiplication by unity. Thus Eq. (C12) can be rewritten 
as 



G C i as5ica i(i) = (+|r e ^ H " (t ' )<it 'oi 

xTe-ti HFp{t ' )dt '0 2 \G). 



(C13) 



With the aid of Eq. ( C5 ) , we can express this quantity 



in terms of the quantum Hamiltonian H instead of Hpp, 

GdassicalW = He^^'Wo W 0l 

xTe-fo H (t')dt' e v N /(2T cl ) ^ G y (C14) 
Rather explicitly multiplying and dividing by y/Zjy, 

e -V N /(2T ci ) 



G 



classica 



i(*) = <+|- 



7- e /o H ^') dt 'OiTe- & H{t ' )dt ' 
x 2 ^/z N e VH ^ 2T ^\G) <C15) 



As discussed in Appendix |A| (in particular, Eq. ( |A20[ )), 
the ground state of the quantum system is given by |0) — 



y/Z^e Vti /( 2T ^\G ). Fu rther invoking Eqs. \C2\CA\ 
can rewrite Eq. (C15) as 



Gclassical(t) = <O|Te/o^(t')rft' 01 

xTe-ti H ^ dt 'O 2 \0). 



(C16) 



Note that in this equation, |0) is the ground state of the 
system defined by the quantum Hamiltonian H . Our re- 
sults above are general. We will shortly use Eq. (C16 1 in 



order to relate it to correlations in the quantum system. 
Under the exchange of t by (it) , the reader may recognize 
Eq. (C16) as a correlation function in the quantum sys- 
tem. One very simple point which is worth emphasizing 
is that not only the ground state of Eq. ( A20 ) is triv- 



ially normalized but, of course, any state formed by the 
evolution with the unitary time dependent time ordered 
exponential U(t) = e~ l -fo #(*')<**'. 

In the quantum arena, it is clear that for a system ini- 
tially prepared in the ground state |0) and then evolved 
with some Hamiltonian H(t) 1 the corresponding correla- 
tion function is given by 

G Quantum (i) = (0|Te , J"o H V )dt ' 1 Te-*/<W* / 

xO 2 |0). (C17) 

By comparing Eq. (C16) with Eq. ( |C17 ), it is seen that 

uantum (t) = G c ,assical(«t). (C18) 

This fundamental identity establishes the connection be- 
tween the overdamped Langevin equation of a classical 
particle at finite temperature and the Schrodinger equa- 
tion of the dual quantum Hamiltonian. We suspect that 
related to this result is the fluctuation-dissipation theo- 
rem that relates correlation functions with the expecta- 
tion value of time-ordered products in equilibrium, see 
for example Ref. IT5l (Chap. 13). 

We next turn to the quantum response function 
^Qu a ntum that monitors the change in the average value 
of 0i as a result of a perturbation 2 . Specifically, we 
consider the Hamiltonian 



H + H 



(C19) 



where H' is a small perturbation which can be expressed 
as H' = — X0 2 - We next review the standard protocol 
for computing the lowest order deviation 



S(O 1 (t)) = (O 1 (t)) x -(O 1 (t)} , 



(C20) 



which we will evaluate within the ground state |0). This 
deviation is readily computed within the interaction pic- 
ture where we evolve with the time ordered exponential 
Texp{-iH't), 

(0 1 (t)) A = ((l-i [ dt! A(i) 2 (t')] O^t) 



x + i 



dt' A(t) 2 {t!) ).(C21) 
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Collecting terms to lowest order, 

a<0i(i)>«i j dt' a(0 (pi(t),o 2 (0]) 

/>00 

= </ dr A(i-r) <[Oi(r),O 3 (0)]) 
Jo 

/oo 
C?T A(*-T) i?Quantum(r). (C22) 
-OO 

As Oi{t) = Te-tti H{t')dt' 0iTe iSlB{t')dt' ^ from the lagt 
line of Eq. (C22|, the quantum response function 



i?Quantum(t) = 10 (t) (0| [Te* # * (*')*' Q x 

xTe- l ^ H ( t ') dt \o 2 ]\0). (C23) 



Comparing Eqs. (C16 C23), we derive Eq.([7]) by further 
expanding it to get the imaginary part of the analytically 
continued classical correlation function, 

-ffQuantum(i) = (G c | assica | (it) - G*| assica | (it)) 

= -2e(<)3G cla55ical ( i <). (C24) 



Appendix D: Arbitrary Number of Time-Ordered 
Operators 

A derivation similar to that in Appendix [C] can be per- 
formed for a correlation function involving an arbitrary 
number of operators. In the classical case, the correlation 
function takes the form of 



G cl = (0 1 (t 1 )O 2 (t 2 )..A,(*»)>, 



(Dl) 



where Oi are arbitrary operators and t\ <t 2 < ... <t n . 

Similar to our calculations in Appendix [C] by Bayes' 
theorem this correlation function is given by 



d dN Xl d dN X 2 ...d dN X n O n 

: ({x n }\Te- HFF{t ' )dt ' \{x n ^})O n ^ 

e - / 9C/ N ({S 1 }) 



.({x2)\Te- ^ Hpp ^ dt ' \{x x }) O x 



J N 



Invoking identity matrix insertions and integrations 
over complete eigensets of states as before, this reduces 
to 

<+|Te-/£ WWOnTe-R^ HFp(t ' )dt 'o„_i 

...re- f *t HFp( - t ' )dt 'Ox\G).(D2) 

Transforming to the quantum Hamiltonian H(t) and 
its respective ground state at time t = yields 

(OlTe-Jt 1 ^Oje"^ H(t ' )dt 'o„_ 1 

...Te-J"*! ff ^*'O!|0). (D3) 



The remainder of the derivation is similar to that in Ap- 
pendix[C] In order to go from the classical to the quantum 
system, we replace the times {£ a }™ = i by {i£ a }„ =1 to go 
from the classical to the many body quantum problem. 



Appendix E: Analytic continuation of classical 
stretched exponentials 

In Appendices [C] and \D\ we explicitly illustrated how a 
Wick type rotation t — » it relates time dependent corre- 
lation functions in the viscous classical systems to those 
in the dual quantum many body theories. For the sake 
of clarity, we explicitly discuss how the analytic continu- 
ation of the classical correlation function should be per- 
formed in some simple yet empirically relevant cases for 
dynamical response functions in viscous classical systems 
wherein Eq. ^ describes the dynamical response. In 
Eq. ([9]), we provided the quantum dual to exponentially 
stretched classical dynamics. In this brief appendix, we 
describe how this result is derived and outline how ana- 
lytic continuations for other response functions of classi- 
cal viscous systems may be analytically continued follow- 
ing such a Wick type rotation. We thank Carl Bender 
for a quick tutorial on aspects of Stokes wedges on which 
this brief appendix heavily relies. 

If G c i assiC ai (t) = J2n B n e^^ T " then, trivially, the 
quantum response function will be uniquely defined and 
given by Rq = B„ cos(i/r„). The same applies, of 
course, for distributions of modes (whence the discrete 
sum over overdamped classical modes n is replaced by 
an integral with some density of modes /(t)). In the 
limit of an infinite number of modes, 



G 



classica 



(*) 



dr7(r')exp(-i/r'), (El) 



with P a distribution that generally is no longer a sum 
of Dirac delta functions. As we will further elaborate 
on in the appendices [see, in particular, Appendix [C] , 

Gquantum(t) G classical(^t) . TllUS, 
p 00 

G Quant u m (t)= dT , f(T , )exp(-it/T / ). (E2) 
Jo 

In the complex r' plane, for any "well-behaved" function 
f(r') that is localized in a region of positive finite r', the 
integral of Eq. (Ell may be performed along any con- 
tour connecting the origin and the r' = 00 along the real 
line such that the contour lies exclusively in the right 
half complex plane (the pertinent Stokes wedge in this 
case) of a positive real component of r', i.e., 5R{r'} > 0. 
Of particular interest to us is the stretched exponential 
form given by G c | ass j ca |(t) = A exp[— (t / r) a ] . Now, in per- 
forming the substitution t — > it to implement the trans- 
formation from Eq. ([El]) to Eq. (E2), we perform the 
rotation 



t = te 



ip : 



Y final 1 



(E3) 
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with if fi na i — 7r(4n+l)/2 where n is an integer. The inte- 
gral of Eq. ( E2 ) remains well-defined in the top complex 



half-plane of a positive real part of i J , i.e., 3?{t'} > (the 
rotated counterpart of the original Stokes wedge) . If ip is 
varied continuously from to 7r/2, there remain contours 
from t' = to r' = oo that appear in the original Stokes 
wedge of !R{t'} > that pass exclusively through the re- 
gion 3{r'} > 0; the integrals along these contours can be 
analytically continued when ip is continuously increased 
from to 7r/2. Thus, we may perform the rotation of Eq. 
( E3 ) continuously increasing ip to represent i as e l7r / 2 and 
replace t — > ie 47r / 2 in the argument of G c i ass i ca i(t). This 
is what we have done in Eq. Q. For other choices of 
n for ip final = 7r(4n + l)/2 as we continuously vary ip 
from its initial value of zero, there will always appear 
situations where the original Stokes wedge will have no 
overlap with its rotated counterpart. Thus, the substitu- 
tion of t — ¥ £e l7r / 2 in the argument of G c i aS sicai (t) is the 
only one that may be implemented out of the possible 
choices in Eq. ( E3 1 in order to evaluate the integral of 



Eq. ( E2 ) . This forms the correct analytic continuation of 
the original real time correlation function of G c i ass i ca i (t) 



the Laplacian is given by V 2 = + ^7^^: it is readily 
verified, as it must self-consistently be, that the corre- 
sponding quantum potential VQ ua ntum = in the region 
outside the range of the interaction. 



2. Harmonic oscillator systems 

As seen by Eq. ([5]), classical systems with harmonic 
potentials Vn map onto quantum systems with similar 
(up to innocuous shifts) harmonic potentials VQuantum- 
That this must be so is readily seen as the ground state 
°f simple quantum harmonic potentials is given by 
a Gaussian. Using Eq. (A20), we see that this indeed 



ofEq. (Ell 



relates to a harmonic classical potential Vn as it must. 
As can be further seen from Eq. (B2), in the case of 
harmonic classical systems, the operators A and A' are 
trivially related to the raising and lowering operators in 
the quantum harmonic problem (and indeed the Gaus- 
sian form of the ground state can, as is very well known, 
be seen from the requirement that the annihilation oper- 
ator must yield zero when acting on the ground state). 



Appendix F: Simple examples of classical to 
quantum correspondence and some aspects 

To elucidate some aspects of the known mapping be- 
tween classical dissipative and quantum systems, we dis- 
cuss several extremely simple examples in d spatial di- 
mensions. 



1. Non-interacting particles 

In the case of a system with a zero potential every- 
where, the quantum ground state wave-function is a con- 
stant in real space. That this is so is seen by our mapping 
as, by Eq.(A20), the classical potential energy vanishes 
everywhere. By Eq. (J5|, the same also occurs for the 
quantum potential which is everywhere zero. 



3. Localization and derealization 

A general localization in the quantum problem will ap- 
pear whenever a corresponding spatial localization occurs 
in the classical system. That is, if the probability of find- 
ing a classical particle outside a distance R from its initial 
location, e.g., the origin, 



h(R) 



d d xP eqml {x) 



(F3) 



x>R 



tends to zero, i.e., h(R) ~ exp(— R/Rq) for R ^> Ro, 
then a localized or bound state is also present in the 
quantum dual. Here x — \x\. In other word s, the part icle 
is localized within a range R. From Eqs. ( A13|A20 |, in 
a system with pair interactions, a localized ground state 
may appear only if 



dvP air (r) 
lim — ^ > 0. 

r— >oo ar 



(F4) 



a. Zero energy bound state 

For a short range attractive potential the zero energy 
eigenstate outside the potential, up to volume normaliza- 
tion factors, given by 



*o(£) 



.4 



1d-2 ' 



Invoking Eq. (A20l, we see that, in this case, 



V^ e ({x}) = 2T d (d-2)ln\x\. 



(Fl) 



(F2) 



Indeed substituting Eq. (F2) into Eq. ([5| and recalling 
that, in its scalar "S-wave" (or u t = 0") representation, 



4. Scaling invariance of time and space 

As is well known, for a homogeneous classical poten- 
tial Vn({x}) which scales as a power (say, p) of the spa- 
tial coordinates \x\, the equations of motion are invariant 
under a simultaneous rescaling of the time coordinates. 
This analysis is typically done for inertial systems. When 
replicated for the over damped system of Eq. (|3| , we find 
that 



Xi — > axi, t — > bt, 



(F5) 



where b plays the role of A and a plays the role of A 1 / 2 
from before, leads to an invariance of Eq. Q if b — a 2 ~ p . 
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By contrast, in the corresponding quantum problem of 
the Schrodinger equation with the Hamiltonian of Eq. 
([5| , a scaling such as that of Eq. (|F5|) is possible only for 



a single case: that of a potential Vn {{x}) that is a log- 
arithmic function of its arguments (or a constant). For 
this particular case, we find that b = a 2 . Correspond- 
ingly, akin to subsection |F 1[ for this particular case, the 
time scales as t ~ 
quantum problem. 



as in diffusion or the free particle 



5. Relation between the classical and quantum 
potentials in the eikonal approximation to the 
Schrodinger equation 

Below we briefly review the eikonal approximation and 
then discuss its relation to the connection between the 
classical and quantum many body potentials as seen in 
Eqs. (5 A20). Towards this end, we write the wavefunc- 



tion as a function of only the phase, the eikonal approx- 
imation [as throughout, we set h= 1], 



* = Ae 



iS 



(F6) 



and substitute this into the Schrodinger equation with 
the Hamiltonian in the second line of Eq. ([5| then we 
will arrive at 



^ 5Z(ViS) 2 + VQuantum(M) + 



as 
at 



2m 



(F7) 



For time independent solutions, ^ = and Eq. (F7) 
rather trivially becomes 

v Quantum ({ a -}) = ^[^v^ 



2m 



(V;S) 2 ]. (F8) 



If we now invoke the correspondence iS <H- — /3Vn/2 then 
Eq. ( F6 1 will transform into Eq. (A20) and , similarly, 
Eq. ( F8 ) will become Eq. ^ relating the quantum po- 
tential energy VQ ua ntum to the classical potential energy 
V N . 



The function \ can take any antisymmetric form. For 
simplicity, we choose it to be a Slater determinant of the 
form 



1 1 



e ik 2 -r 1 e ik 2 -r 2 



gifcjv-fi e ik N -f 2 



ik 2 -r N 



a ik N -r N 



, (G2) 



with Q the volume of the system. 
As idt^o = H^ , 



id t (^ox) = x(*$*o) + *o(*$x) 

= X H^0 + * ^SIaterX, 



(G3) 

(G4) 



where -Esiater is the energy of the free particle system 
described by %. 



uantum 



(G5) 



The potential energy operator leads to xC^Quantum^o)- 
The kinetic energy operator T generates three terms of, 
respectively, the forms Ea( V aX)'0 O , Ea( v aX) ' (y a tp°), 
and X)a x(VaV J )- The first and the last terms represent 
the term proportional to -Bsiater and the original bosonic 
kinetic energy respectively. The second term, that of the 
mixed gradients, is proportional to ~^2 a k a . For a sys- 
tem invariant under parity, this sum vanishes. Up to 
an innocuous phase factor, the evolution given an ini- 
tial fermionic wavefunction of Eq. (Gil, will be thus 



identical to that with the bosonic wavefunction and 
all correlation functions will be identical to those which 
we earlier computed for the bosonic system. That is, 
the general time dependent correlation functions given 
an initial fermionic state of Eqs.( A20|G1 ) will adhere to 
the general t —} it rule which we detailed in earlier sec- 
tions. A notable difference by comparison to the bosonic 
case, however, is that the wavefunction of Eq. (Gl) at 



an initial (or at a final) time is, generally, not a ground 
state of the Hamiltonian H. 



Appendix G: Slater-Jastrow forms 

The general results presented thus far may, in some 
instances, be generalized to describe fermions. A lim- 
ited extent extension is the one concerning the evolution 
starting off from an initial Slater-Jastrow type fermionic 
wavefunction. As we have emphasized earlier, if the 



in Eq. (AlOl are symmetric under the exchange of i and 
j, the resulting wavefunction obeys Bose statistics. This 
symmetry is maintained for the ground state as it is a 
Jastrow function given by Eq. (A20). Fermionic wave- 



functions are afforded by the product of the symmetric 
boson ground state and an antisymmetric term, 



(Gl) 



Appendix H: Extensions: complex wavefunctions 
and lattice systems 

We now briefly suggest and elaborate on several exten- 
sions of our calculations thus far. We will illustrate and 
suggest how our results may hold for (i) several lattice 
extensions as well as (ii) general systems with complex 
wavefunctions. This latter case allows us to go from a 
given quantum mechanical problem to a corresponding 
classical one. 

(i) Lattice systems. If we replace Eq. |3]) by 



7i 



S 



dt 



Vn(0i, 0jv) + rji(t), 



(HI) 
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to describe a classical lattice system with fields <fii at the 
various lattice points i, then trivially replicating all of our 
calculations thus far with the exchange Xj — > 4>i (includ- 
ing in all gradient or variational derivative operators), we 
will arrive at a corresponding quantum lattice system. 

(ii) Complex wavefunctions. The similarity transfor- 



mation of Eqs. (A9 A19 A21 1 captures a simple mathe- 



matical identity between the generalized probability dis- 
tribution of a classical system, obeying the Fok ker-Planck 
dynamics with an operator Hpp (Eq. (A7|), and the 
wavefunction obeying the Schrodinger equation of the 
quantum dual Hamiltonian H . Given this relation, it 
is possible to, formally, consider extensions in which the 
function evolving with the Fokker-Planck dynamics need 
not be a probability distribution. Most of our results 
concerning temporal correlations may hold under such a 
generalized interchange if we invoke Eq. (A20) to de- 



fine (when given a quantum problem) a corresponding 
classical system which need not be a bona fide physi- 
cal Boltzmann distribution as in the scalar bosonic sys- 



tems which we primarily focused on thus far and employ 
| "Jo 1 2 as the initial (or final) time weight in the multi- 
ple time classical correlation function of Appendix [C] In 
appendix |F 5| we examined a formally imaginary coun- 
terpart to Vjvj and explicitly demonstrated how it leads 
to standard results. Thus, given a quantum wavefunc- 
tion, we may consider its logarithm to correspond to a 
classical potential Vn. Wavefunctions of spinless Fermi 
systems cannot be purely positive and for these com- 
plex (as well as divergent) potentials will formally arise. 
There may be subtleties however in our imaginary time 
(t — > it) analytic continuations when Vn is not purely 
real (and the system effectively not purely dissipative) 
which are more complex than those which we invoked 
thus far in our analysis of real Vn which led to response 
functions of pure damped modes and their superpositions 
such as those which we encountered in Eq. ^ [see also 
Appendix [e] . Physically, these are related to analogs of 
classical systems with instantons and tunneling events 
(the behavior for the pure dissipative system) when these 
further exhibit non-damped oscillatory behavior. 
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